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In a recent experiment [PRL 100, 040404 (2008)] an analog of photon-assisted tunneling has been 
observed for a Bose-Einstein condensate in an optical lattice subject to a constant force plus a 
sinusoidal shaking. Contrary to previous theoretical predictions, the width of the condensate was 
measured to be proportional to the square of the effective tunneling matrix element, rather than 
a linear dependence. For a simple model of two interacting bosons in a one-dimensional optical 
lattice, both analytical and numerical calculations indicate that such a transition from a linear to 
a quadratic dependence can be interpreted through the ballistic transport and the corresponding 
exact dispersion relation of bound boson pairs. 
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I. INTRODUCTION 

Bose-Einstein condensates (BECs) in optical lattices 
provide an excellent tool to study solid state systems JJ] ■ 
One of the methods which are currently established ex- 
perimentally d, H, 13] for BECs in an optical lattice 
is tunneling control via time-periodic potential differ- 
ences [1, 0: S IB • Effects investigated theoretically in pe- 
riodicall y sh aken systems include multi-particle entangle- 
ment [9|, [10| and nonlinear Landau- Zener processes 

The experimental realization 0, [l^, EB of destruction 
of tunneling via time-periodic potential differences @ 
was the breakthrough for tunneling control via time- 
periodic potential differences. The systems used so far 
in experiments are as divers as BECs in an optical lat- 
tice m, single particles in a double well and light in 
a double- well system fTs"] . 

Reference [8] suggested to use time-periodic shaking 
for a tilted double-well potential to measure an analog 
of photon- assisted tunneling - the "photon" -frequencies 
corresponding to shaking-frequencies in the kilo-Hertz 
regime. Recently, photon-assisted tunneling was ob- 
served for a BEC in a periodically shaken optical lat- 
tice Q. When the BEC was allowed to expand for some 
time, the width of the condensate was measured to be 
roughly proportional to the tunneling matrix element for 
small BECs whereas it was proportional to the square of 
the tunneling matrix element for larger condensates. 

While this might be an indication of a transition from 
ballistic to diffusive transport Q , the verification of such 
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a transition requires the precise measurement of the time- 
dependence of the width of the wave function in order to 
distinguish both regimes of transport (see, e.g., Ref. [T^ 
and references therein). Without such measurements, 
other explanations for the observed dependence of the 
width of the BEC cannot be excluded. 

Here we develop an alternative interpretation of the 
experimental results: Using two interacting bosons in a 
simple tight-binding one-band model we show that the 
transition from a linear to a quadratic dependence on 
the tunneling amplitude might be an interaction-induced 
effect within the ballistic regime. We demonstrate, both 
with the help of analytical arguments and numerical sim- 
ulations, that the wave function of an interacting boson 
pair indeed behaves qualitatively very similar to the BEC 
in the experiment: On the one hand, weakly interacting 
particles can reproduce the linear scaling with the tunnel- 
ing matrix element observed for small condensates. On 
the other hand, more strongly interacting particles repro- 
duce the quadratic scaling for larger condensates (which 
have a larger total interaction energy). 

The paper is organized as follows. In Sec. |TT] we in- 
troduce the one-dimensional Bose-Hubbard Hamiltonian 
with a constant force and time-periodic driving used to 
model the experimental situation. In the case of high 
driving frequencies, the time-dependent Hamiltonian can 
be replaced by an effective, time-independent Hamilto- 
nian. For this Hamiltonian, exact two-particle energy 
eigenstates describing bound boson pairs can be derived 
[ill [l6| as is explained in Sec. IIIII In Sec. II VI we investi- 
gate the dynamical behavior of two-particle wave pack- 
ets. Employing the obtained exact dispersion relation 
for interacting boson pairs, we derive the dependence of 
the width of the wave packet after some time of free ex- 
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pansion on the ratio of interaction and tunneling matrix 
element. These theoretical considerations are supported 
by numerical simulations of the two-particle Schrodinger 
equation for the full driven Bose-Hubbard Hamiltonian. 
Some conclusions are drawn in Sec. fVl 



Without interactions (17 = 0), the effective Hamilto- 
nian reduces to a well-known single-particle Hamiltonian 
with extended Bloch-waves as eigenfunctions 



Jkjdi_ 



|J>, 



(5) 



II. THE MODEL 

Using the notation of Ref. the one-dimensional 
Bose-Hubbard Hamiltonian with time-periodic shaking 
describing the experiment can be written as: 



Ho = --/Ef^S+i 



-AE E jfi-j + K cos{ujt) E jfij J 



1) 



(1) 



where the operators c'p annihilate (create) bosons at the 

lattice site j and fij = cjc^ are number operators; J is 
the hopping matrix element, AE = FcLl is the potential 
difference of two adjacent wells with lattice spacing cJl, 
w / (27r) the frequency with which the system is shaken 
and K the amplitude of the shaking. Thus, the Hamil- 
tonian ll]) describes a system of interacting bosons in 
a tilted and driven optical lattice. For an untilted and 
undriven lattice we have AE — and K — 0. 

Floquet-theory can be applied to understand the 
physics of such a driven system. For not too low driving 
frequencies (hv ^ J), the resonance condition 



nhuj = AE 



(2) 



with integer n leads to photon-assisted tunneling ([1,01). 
In the high-frequency limit (Sw ^ J and fuj ^ U |18| l 
many aspects of the physics behind the system can be 
understood by replacing the time-dependent Hamiltonian 
with constant force by a time-independent Hamiltonian 
without any additional force: 

^off = - Jeff E + + E " ^) ' 

3 j 

(3) 

where the effective tunneling matrix element is given by 



_ K 
njjj 



(4) 



with the nth order Bessel function J^n- Thus, in this 
time-periodic system photon assisted tunneling can lead 
to an increase of the tunneling amplitude (as the energy 
difference between neighboring wells is removed in the 
effective Hamiltonian). Moreover, within an n-photon 
resonance, tunneling can also be suppressed by tuning 
the driving amplitude K such that the ratio Kq of driving 
amplitude and huj corresponds to a zero of the Bessel 
function J^n{Ko). 



where \j) is the Wannier-function at lattice-site j, with 
the corresponding energy eigenvalues 



E{k) — — 2Joff cos(fco?L)- 



(6) 



Within a parameter regime for which the cos-dispersion 
relation can be replaced by a quadratic dispersion rela- 
tion, E{k) ~ — 2Joff + Jcsk'^dl, the single particle in a 
tight-binding lattice behaves like a free particle. In fact, 
such model Hamiltonians can be used to do numerics for 
a free particle by setting 



Jcffd? = 



2m ' 



(7) 



where m is the mass of the free particle with dispersion 
relation Ef{k) = h^k'^/{2m). 

Without interaction, a BEC of N bosons would simply 
be the product of single-particle wave functions. If 
one measures the width of a condensate after a certain 
time of free expansion Q, this thus corresponds to the 
popular text-book exercise [1^] of calculating the width 
Ax{t) for a single free particle. Starting from a Gaussian 
wave packet at time zero. 



ip{x,t = 0) 



1 



(27ra2)i/4 



exp 



X 



finds 



Ax{t) 




ht 



(8) 



(9) 



where Ax{tY = {tP{t)\x'^\'iP{t)) ~ {tP{t)\x\il){t))'^ . 

For a non-interacting BEC in an optical lattice one 
can thus expect to find for not too small free expansion 
times t: 



Ax{t) cx \E"{Q)\t, 



(10) 



where the dashes denote derivatives with respect to the 
argument (fc). This equation was used in Ref. Q to mea- 
sure the effective tunneling matrix element Jcs as using 
Eqs. and HO]) one has: 



Ax{t) cx 2| JoffML*- 



(11) 



In order to see if this relation always survives interaction, 
two interacting bosons are investigated in the following 
sections, starting with the construction of correspond- 
ing energy eigenstates. Two-particle effects in optical 
lattices are interesting both experimentally and theoret- 
ically [15|. Recent investigations include superexchange 
interactions for mixtures of different spins [2^. Moti- 
vated by the experiment ^ , this manuscript concentrates 
on pairs of indistinguishable particles. 
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III. EXACT ANALYTIC TWO-PARTICLE 
EIGENFUNCTIONS 



IV. BALLISTIC EXPANSION OF 
TWO-PARTICLE WAVE PACKETS 



The aim is to find exact analytical expressions for the 
eigenf unctions of the Haniiltonian ([3]) for two bosons and 
non-zero interaction [/ 7^ 0. Rather than using the ap- 
proach via Green's functions of Ref. 15], one can proceed 
along the lines of Ref. [21] to show in a straightforward 
calculation (see the appendix for details) that a large 
class of two-particle wave functions is given by 



E 



a^.,^,{k)\v)\^l), 



(12) 



where v < required because the bosons are indistin- 
guishable, and \v)\v) corresponds to the Fock state with 
two particles at lattice site v. The coefficients a^^^i are 
given by 



{k) 



where 



^^(fc) = (r/x-)!'" "I exp {ikd\,{v + 



-1 
+1 



[//Jeff > 0, 

U/J,ff < 0, 



and 



C/2 



\u\ 



16 J^g cos'^lkdh) 



|4Jcff| cos(fcdL) 



(13) 



(14) 



(15) 



(16) 



As for the Bloch-waves ([S]) and for plane waves for free 
particles, these wave functions cannot be normalized in 
the usual sense. Nevertheless, to avoid divergence for 

— iy\ ^ 00, one needs |a;_| < 1 and thus cos(fcdL) > 0. 

The energy eigenvalue of the state \(f>k) can be written 
as (see Eq. and Ref. 0): 



E2{k) = -4?7Joff. 



[/2 



+ cos2(Ml), (17) 



and one thus obtains 



16 Je^ffdE 



(18) 



The solutions calculated above represent the bound 
states of a boson pair. In addition to these solutions one 
has a continuum of scattering eigenstates [12, which 
will play, however, only a minor role in the description of 
the physical situation considered here: The initial wave 
function is given by a narrow Gaussian centered around 
a lattice site. Thus, both particles are likely to be sitting 
at the same lattice site and, therefore, the initial state 
lies nearly entirely in the subspace of bound states of the 
two-boson system. 



The value that is relevant for the spreading of two- 
particle wave packets is not the energy E2{k) of a pair but 
the energy per particle. Hence, for the effective Hamil- 
tonian Q we have 



\E"m 



(19) 



Equations (|T0|) and (|T9|) therefore show that the width of 
the BEG in the driven and tilted lattice is given by 



Aa; j^ff {t) cx 



Vl64^ + ^' 
from which we find the following limiting behavior. 

2\.hs\dlt : |C//Jcff| «1, 



Axj^,, (i) (X 



^Jls/\U\dlt : \UIM » 1, 



(20) 



(21) 



where Jcg = JJn{Ko) was introduced in Eq. The 
corresponding expression for the width in the undriven 
and untilted lattice, which we denote by Axj(i), is ob- 
tained by replacing Jcb by J, 



Axj{t) 



SJ^dLt 

V16J2 + C/2' 



(22) 



Thus we see that the transition from a linear to a 
quadratic dependence on the Bessel function observed 
in Ref. 3] can, within this simple two-particle model, be 
explained as being a continuous transition based on the 
dispersion relation (|17p for two interacting particles. 

It is important to note that one does not have to wait 
until the dependence of Ax(t) on t becomes linear in 
order to see the scaling when comparing, e.g., an un- 
driven system without tilt with a periodically shaken, 
tilted system with n-photon-assisted tunneling. In the 
high-frequency limit, the main difference is that the mod- 
ulus of Jeff will be lower than J by a factor of | Jn {Kq ) | . If 
one plots the width of an initially localized wave packet 
(cf. Ref. 3) as a function of r = Jt/h, the undriven 
system will thus spread faster. Only by rescaling the 
time scale for the undriven system one can hope to make 
both functions agree. For weak interactions, the follow- 
ing equation relates the width of the undriven system at 
time \Jn{Ko)\T to that of the driven system at time r: 



AXJ {\Jn{Ko)\T) ^ Axj^^^{t) , 



u 

JcS 



< 1. (23) 



For stronger interactions another factor of \JniKo)\ is 
necessary: 



Axj {\JJKo)\^t) ~ Axj^,,{t) , 



U 



Jeff 



> 1. (24) 



The above reasoning explains why the width of the 
BEG can be proportional to J^{Kq) even within the 
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ballistic regime. However, the experimental results of 
Ref. [31] show even more, namely that Axj^j,. (t) is not 
only proportional to J^{Kq), but also that the ratio of 
the width in the tilted and driven lattice to the width in 
the untilted and undriven lattice is approximately equal 
to J^iKo), i. e. 

Axj^,, {t) 2 

-^-^^J^iKo). (25) 

We demonstrate in Figs. [1] and [2] that it is indeed possible 
to find parameters which reproduce the experimentally 
observed behavior within the framework of the theory 
presented here. In Fig. [T] we compare the driven inter- 
acting system {U ^ 0) with the undriven system for non- 
interacting particles {U = 0). We see that for U j J ^ 3.6 
(dotted curves) the left-hand side of Eq. (|25|) nearly lies 
on top of the curves representing the right-hand side of 
this equation (dashed curves). For lower interactions 
(dash-dotted curves), the scahng is again cx | Jn(-ft'o)| as 
in the single-particle case. 

In Fig. [1] the width of the wave function was compared 
to the width for non- interacting particles. However, the 
observed scaling even occurs (for larger interactions than 
those chosen in Fig. [T]) when comparing the untilted, un- 
driven interacting system with the periodically driven in- 
teracting system: Figure [2] shows that again a very sim- 
ilar scaling of the width of the condensate is found. 

Although several papers have shown the validity of 
the effective Hamiltonian approach used so far (see, e.g., 
Ref. [l^ and references therein), one should demon- 
strate that it also is valid for the present situation. 
To do this, we numerically solve the time-dependent 
Schrodinger equation corresponding to the full Hamilto- 
nian ([1]) (Fig. [3]) . As the initial wave function we choose 
a Gaussian function for the ccnter-of-mass wave function 
(corresponding to an initial confinement via a harmonic 
trapping potential), 

|V'(i = 0)) =^exp(-a2(fc_fc^)2/2) jv^^), (26) 
fe 

where the sum over all possible fc values (rather than 
an integral) is necessary as numeric calculations can- 
not be done in infinite lattices. The initial wave func- 
tion is sitting in the middle of the lattice (for an odd 
number A^l of lattice-sites, the sites can be labelled as 
j = -(A^L - l)/2 . . . (A^L ~ l)/2). In the finite, shaken 
lattice relevant for the numerics in this paper, vanishing 
boundary conditions are a suitable choice [2j. For the 
center-of-mass part of the wave function (cf. Eq. (|13p ). 
Mc.o.m. = (/^ + '^)/2, the wave- vector is 2fc, possible val- 
ues for k are thus n7r/[(A^L ~ l)c?L]7?^ = 1,2,3, . . .. The 
initial wave function was calculated without any initial 
momentum (fco — 0) and with a — lOdL- For the un- 
driven system, U / J = 10.0 was chosen in the initial wave 
function, for the driven system with Kq — 2.0, we chose 
U I J = 10.0/j7i(2.0) to mimic an experimental situation 
where the initial wave function is prepared in a harmonic- 
oscillator potential and the periodic shaking is switched 
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FIG. 1: The width of the two-particle wave function as a 
function of the ratio of driving amplitude and driving fre- 
quency A'o = K/ifiULi). The width is normalized to the value 
for non-interacting particles {U — 0) in an undriven, un- 
tilted lattice. The curves are based on the analytic Eqs. Q, 
(|20p and H22p . For the one-photon resonance (upper panel) 
one has Jeg = JJi{Kq) and for the two-photon resonance 
(lower panel) Jcs ~ JJ2{Ko)- Solid curves: |^7n(ifo)|, dashed 
curves: J'^{Ko), dash-dotted curves: U/J = 0.6, dotted 
curves: U/J = 3.6. 

on before switching off this potential. The quadratic scal- 
ing as predicted in Eq. ([M]) can thus indeed be observed 
numerically (Fig. [3]). 

V. CONCLUSION 

Photon-assisted tunneling in a periodically shaken op- 
tical lattice was investigated for two interacting bosons. 
Both numerical and analytical calculations were done 
for a periodically driven one-band Bose-Hubbard model. 
Figures demonstrate that the experimentally ob- 
served [3] dependence of the width of the BEG on the 
square of the tunneling matrix element can be explained 
- at least qualitatively within the simple two-particle 
model investigated here - as being an interaction-induced 
effect which is based on the dispersion relation for bound 
boson pairs. 

While our simplified approach thus explains some as- 
pects of the experiment (3| , calculations for larger parti- 
cle numbers (which are not straightforward to generalize 
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FIG. 2: The width of the two-particle wave function as a func- 
tion of the ratio of driving amplitude and driving frequency 
Ko = K/{hLu) normalized to the value for interacting parti- 
cles in an undriven, untilted lattice. Labels can be found in 
Fig. [1] but for the dotted curves which are calculated for an 
interaction of U/ J = 10.0. As in Fig. [T] the transition from 
linear to quadratic dependence on Jn{Ko), n = 1,2, can be 
observed. 



from the method presented here) are likely to lead to 
further insights into the experiment for which additional 
effects like decoherence by particle losses might also play 
an important role. An experimental measurement of the 
time-dependence of the width of the wave function would 
be of great interest for the theoretical analysis and the 
modelling of the transport properties of quantum con- 
densates. In particular, it would be interesting to see if 
indeed a transition from ballistic to diffusive transport 
takes place, or if the simple two-particle model presented 
in this paper can explain the relevant features of the ex- 
periment. 
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FIG. 3: Numeric simulation of a two-atom wave function in 
an optical lattice with TVl = 201 lattice sites in the high- 
frequency limit (huj — 40 J) with large interactions U/J — 
10.0 (cf. Fig. [2)) and initial wave functions given by Eq. (I26|l . 
a) The width is plotted as a function of dimensionless time 
r = Jt/h. Dotted curves: the undriven, untilted lattice. 
Dashed curves: the same lattice but with periodic driving for 
parameters corresponding to a one-photon resonance @] 
{Ko — 2, cf. Eq. Q). The wave packet in the undriven 
system spreads faster than in the driven system, b) With- 
out the knowledge of Figs. [1] and (2] one might expect both 
curves to agree if one plots Aa;(ri), ri = \J'i{2)\t for the un- 
driven system. However, one needs another factor of Ji{2): 
for T2 = [J7i(2)]^r, the data of the undriven system plotted 
as Ax{t2) coincides with Ax{t) for the driven system (c). 



APPENDIX A: EIGENFUNCTIONS 



In order to show that exact two-particle eigenfunc- 
tions [23] of the effective Hamiltonian ([3]) are indeed given 
by Eqs. (fT2|) and ([13]), one can start with 



j 
j 

j 



%/2|m-1)|m), (A1) 
V2\fi)\^i + 1) , 

\^i~2)\^i) + V2\^i-l)\^l-l) , 

|Ai-l)|/i+l)+y2|/i>|Ai) , 



6 



and for v < ii — 1 

J2c]c^^,WM = W-iM + , (A2) 

Using the notation of Eq. (fT2)) one thus has 

Hes\^k) = uY,'^^^Ak)\^^)\^^) (A3) 
/J 

-Jeff ^ x/2a^^^(fc) (l/i - + \^I)\fi + 1)) 
-Jeff a^^^ijk) (Im - 2) \n) + V2|Ai - 1) - 1) 

-Jeff 5Z «A.-l,M(fc) (l-" - 1)1/^ + 1) + 

-Jeff J2 ('•^Ak) iW - + - I)) 

v<ix—l 

-Jeff 51 + + 

In oder to show that this indeed leads to an eigenfunction 
of the effective Hamihonian ([3]), we use Eq. and 
start with the last two lines of Eq. (jA3p . Performing 
appropriate shifts of the summation indices these lines 
can be written as 



- Jeff^(& 

- Jeff J2 iK+iAk) + bu.t.-i{k))W)\tJ.). (A4) 

The second, third and fourth line of Eq. (lASp can be 
combined to yield 

- Jeff I/— 1,1/ 

- Jeff ^ (6t.+i,p(/c) + 

- Jeff ^ (A5) 



Adding (|A4p and (|A5p and including the first line of 
Eq. jMl) we find 

i?eff|0fc) = -^^6,,,(fc)|i.)|i.) (A6) 

- -^Cff 

- Jeff ^ + 6,._i,^(fc)) 

- Jeff ^ (Wi.A'Cfc) + K,iJ,-i{k)) 



In order to simplify the last two lines of Eq. (jA6p we use 
the notation of Eqs. p4p - (fTB|) to obtain 

b„+i,^{k) + 6,.-i,^(fc) + b^^i_,+i{k) + 6i.,p-i(fc) 



bvAk)^ 



b^^^{k)2r] cos (kdi,) [x- +x_^] 



with 



[/2 



+ 1 



1/2 



,16j2gC0s2(fcdL) 

Thus, for V < ^ we have 



{llMHefiW) (A7) 
/ 2 \l/2 

= -4??JeffC0s(fcdL) ( i6ja,e^os^(fcdL) + 

The terms with = /i in Eq. (jA6p yield: 

(C//V2)6,,,(fc) - JeffV2[6,^,+i(fc) + (fc)] 
= C//\/2 — Jcff\/2x_2?7Cos(fcdL) bi,,y{k) 
= [C/ - Jcff4a;_77C0s(A:(iL)] ai,,^(A:), 

and thus (cf. Eq. (dH])): 



(zy|(z.|iJeff|0fe) 

= -4?7Jcff cos(/cdL) (y6J^ 



U + 77 Jeff 



J 



eff 



COS^(/crfL) 



1/2 

1) a^Ak) 



(A8) 



where f/ + ?7 Jeff |C/|/| Jeff I = [Eq. US])]. This shows that 
\4)k) indeed is a 2-particle eigenfunction of the effective 
Hamiltonian, 

H,s\ct'k)^E2{k)\4>u), 



with the energy eigenvalue 



E2{k) = -4?7Jeff. 



f/2 



16Je'ff 



cos^ [kdi 



(A9) 



Note that cos(fcdL) was required to be positive for 
to have a modulus lower than or equal to one, because 
otherwise the wavefunction would diverge. 



l'<jLl 
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